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Lattice QCD simulations tend to get stuck in a single topological sector at fine lattice spacing, 
or when using chirally symmetric quarks. In such cases computed observables differ from their 
full QCD counterparts by finite volume corrections, which need to be understood on a quantitative 
level. We extend a known relation from the literature between hadron masses at fixed and at un¬ 
fixed topology by incorporating in addition to topological finite volume effects, also ordinary finite 
volume effects. We present numerical results for SU(2) Yang-Mills theory. 
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1. Introduction 

In QCD simulations at small lattice spacings a<0.05fm, algorithms typically have severe prob¬ 
lems in generating transitions between different topological sectors. This problem of topology 
freezing is expected to appear for any lattice discretization of the quark and gluon fields fflGl. 
For certain discretizations, e.g. chirally symmetric quarks, this problem is even present on coarser 
lattices 001 • In specific cases it might be motivated to fix topology on purpose. For example, 
when using a mixed action setup with light overlap valence and Wilson sea quarks, one observes a 
rather ill-behaved continuum limit 0|0. This is due to near-zero modes of the Dirac operator in 
the valence sector, which are not compensated by corresponding modes in the sea. A possibility to 
circumvent this imbalance could be to restrict the lattice simulation to topological charge Q = 0, 
where such near-zero modes are absent, e.g. by employing topology conserving actions [7-5). 

Methods to extract physically meaningful results from simulations at frozen or fixed topology 
have been proposed iflOl 1JJ and tested [12-22] in various models and theories. In this work we 
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extend these methods by also including ordinary finite volume effects. Such a combined treatment 
of both fixed topology and ordinary finite volume corrections is expected to be particularly impor¬ 
tant for QCD at light u/d quark masses. We test our equations in SU(2) Yang-Mills theory at fixed 
topology. 


2. Topological finite volume effects 

In iflOl 1JJ an equation has been derived relating a hadron mass Mq y obtained at fixed topo¬ 
logical charge Q and finite volume V to its counterpart M at unfixed topology (i.e. the physically 
meaningful hadron mass), 

Mqv =M+—!— ) +e( k- Y (2.1) 

e 2%v V XtVj \itiV) 2 ) 

where M" denotes the second derivative of M with respect to the 6 angle at 9 = 0, and Xt the 
topological susceptibility. This equation illustrates that fixed topology corrections are finite volume 
effects, i.e. effects suppressed by inverse powers of V. It is straightforward to extract a physical 
hadron mass M from computations at fixed topology: one just has to fit eq. (12.11) to the available 
fixed topology and finite volume hadron masses Mqy, where M, M" and Xt are the fit parameters 
(examples of this procedure can be found in ll2j-[2l]|). 

In Fig. |T| we show recent results for SU(2) Yang-Mills theory (standard plaquette action, gauge 
coupling j8 = 2.5, i.e. lattice spacing a « 0.073 fm II23T0 . The static potential ''Vqq&y at separation 
r — 6a (which can be interpreted as a mass), has been computed in different topological sectors with 
topological charges \Q\ = 0,1,2,3,4 and for different volumes V = 14 4 ,15 4 ,16 4 ,18 4 (4000 gauge 
link configurations have been generated for each of the four volumes). Note that the discrepancies 
between static potential results "Vqqny (r = 6a) at different topological charges Q are clearly visible, 
in particular for small volumes vH This demonstrates the necessity of using specific methods to 
properly deal with topological finite volume effects. The curves represent a single fit of eq. (12.11) to 
the lattice static potential results 'Vqq&y (;r = 6a) . The fit is of good quality, i.e. x 2 $ 1 • The resulting 
y q q{r = 6a) — 0.3097(5) is in excellent agreement with a corresponding standard computation at 
unfixed topology, which yields 7^(r = 6a) = 0.3101(3). 

3. Ordinary finite volume effects 

Usually lattice simulations are performed at finite spatial volume L 3 with periodic boundary 
conditions. Consequently, a hadron at x will interact with images of itself, e.g. at x±Le x , xiLe^ 
or x±Le z . Such interactions cause a shift in the hadron mass compared to infinite spatial volume, 
as first derived in lf25l . 

The corresponding equation to describe these ordinary finite volume corrections (i.e. finite 
volume corrections not related to fixed topology) of the static potential in Yang-Mills theory is 

M(L)-M(L4»)ociexp|-^j, (3.1) 


1 Similar observations for the pion mass have been reported in HD. 
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Fig. 1. y q q,Qy(r = 6a) as a function of 1 /V. The curves represent the fit of the lattice static potential results 
to eq. (12,1b . 


where M = y qq {r) and m is the mass of the lightest particle, i.e. the J PC — 0 ++ glueball. In Fig. [2| 
we confront this equation with lattice SU(2) Yang-Mills results for y qq {r = 3 a) and find excellent 
agreement (again we have used j8 = 2.5 and generated 4000 gauge link configurations for each 
of the eight volumes V = 10 4 , ll 4 ,12 4 ,13 4 ,14 4 ,15 4 ,16 4 ,18 4 ). For L > 14 ordinary finite volume 
effects are negligible. For smaller L, however, there are sizeable corrections, which need to be 
taken into account, in particular when using such volumes for computations at fixed topology as 
presented in the previous section. From the fit of eq. (13.1b to the lattice results shown in Fig. [2l one 
can even extract the J PC — 0 ++ glueball mass with remarkable precision, m = 0.74(4). This is in 
perfect agreement with the result obtained by a standard lattice computation of a glueball 2-point 
correlation function, m = 0.723(23) [2611. 

4. Combining topological and ordinary finite volume effects 

In the SU(2) example discussed in Section [2| it has been possible to analyze fixed topology 
results using eq. (12.1b in a meaningful way, i.e. without taking ordinary finite volume effects into 
account. Since the mass of the lightest particle, the J PC = 0 ++ glueball, is quite large, ordinary 
finite volume effects are strongly suppressed for large volumes. As indicated by Fig.[2]and as done 
in Section O one just has to discard volumes with L < 14. In QCD the situation is expected to 
be more difficult, because there the lightest particle, the pion, is much lighter than the J PC = 0 ++ 
glueball of SU(2) Yang-Mills theory. Moreover, lattice simulations of QCD, in particular at large 
volumes, are extremely demanding with respect to high performance computer resources. There¬ 
fore, it is highly desirable to combine eqs. (12.1b and (13.1b . i.e. to obtain an expression describing 
both topological and ordinary finite volume corrections to hadron masses. 

To derive such an expression, one has to consider ordinary finite volume effects also at non- 
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Fig. 2. The dependence of y q q{r = 3a) on the periodic spatial extension L of the lattice (at unfixed topology). 


vanishing 6 angles using equations analogous to (13.11) . These equations are the starting point for 
a lengthy calculation similar to that leading to eq. (12.11) (a detailed derivation of eq. (12.11) can e.g. 
found in lfT8l ). The resulting expression describing both topological and ordinary finite volume 
effects takes the form 


M QV= M +^M , '{\ 

-df.+J 

L V 11,V V A 






(4.1) 


where ordinary finite volume effects for the topological susceptibility % t have been neglected, since 
they are expected to be tiny [27,28]. Note that in addition to the parameters M, M" and % u which 
are already present in eq. (12.11) . there are four more parameters, m, m", A and A", characterizing 
combined topological and ordinary finite volume corrections. 

In Fig. [3] (top) we show a plot similar to that from Fig. [lj this time for y qq ,Qy{r = 3a). More¬ 
over, also results for small volumes V = 11 4 ,12 4 ,13 4 are included. The curves correspond to 
eq. (12.11) with the fit parameters M, M" and Xt determined by a fit to the large volumes V = 
14 4 ,15 4 ,16 4 ,18 4 , where ordinary finite volume effects are negligible. There is a strong discrep¬ 
ancy between these curves and the lattice results for the small volumes V = ll 4 ,12 4 ,13 4 . This is 
expected, since ordinary finite volume corrections are not part of eq. (12.11) in particular for small Q. 



^qq,Qy( r — ^ a ) 

in 

Xt x 10 5 

fit results, eq. (14.11) 

unfixed topology results If26li29l 

0.16437(15) 

0.16455(7) 

0.67(10) 

0.723(23) 

9.5(2.0) 

7.0(0.9) 


Table 1. Results for the static potential y qq (r = 3a), the mass m of the J PC = 0 ++ glueball, and the topolog¬ 
ical susceptibility %, obtained by a fit of eq. (14.11) to fixed topology lattice results yqq,Qy( r — 3^)- 
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Fig. 3. y q q,Qy(r = 3a) as a function of l/V. (top): The curves represent the fit of eq. (12.1b to the lattice 
static potential results for large volumes V = 14 4 ,15 4 ,16 4 ,18 4 . There is a strong discrepancy between these 
curves and the lattice results for the small volumes V = 11 4 ,12 4 ,13 4 . (bottom): The curves represent the 
fit of eq. (14.11) to the lattice static potential results for all volumes V = 11 4 ,..., 18 4 . There is almost perfect 
agreement, even at small volumes and <2 = 0. 


In Fig. [3] (bottom) we show the same lattice results for y qq &y(r = 3a). This time, however, the 
curves correspond to eq. (14.11) with the fit parameters M, M", % u m , m", A and A" determined by a fit 
to all seven volumes V = ll 4 ,12 4 ,13 4 ,14 4 ,15 4 ,16 4 ,18 4 . There is almost perfect agreement, even 
at small volumes and for Q — 0. The extracted “hadron mass” y qq &y(r = 3 a) is consistent with 
a corresponding computation at unfixed topology and also the glueball mass m and the topological 
susceptibility % obtained by the fit are in fair agreement with reference values, cf. Table [U 

To conclude, we have incorporated ordinary finite volume corrections into an existing relation 
between hadron masses at fixed topology and physical hadron masses (i.e. hadron masses at un¬ 
fixed topology). We have successfully tested the resulting equation in SU(2) Yang-Mills theory by 
studying the static potential at fixed topology. As an outlook, we plan to extend these tests to QCD 
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in the near future. 
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